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^ . ABSTRACT 
O 

CN ■ We present the results of a systematic analysis of the world sample of 

Ph 1 optical/near-infrared afterglow light curves observed in the pie-Swift era by the 

■ end of 2004. After selecting the best observed 16 afterglows with well-sampled 

light curves that can be described by a Beuermann equation, we explore the 

parameter space of the light curve parameters and physical quantities related 

CN ■ to them. In addition, we search for correlations between these parameters and 

• the corresponding gamma-ray data, and we use our data set to look for a fine 

ON ■ structure in the light curves. 
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Nearly ten years after the discovery of the first Gamma-Ray Burst (GRB) afterglow in 
the optical/near-infrared (Groot et al. 1997; van Paradijs et al. 1997), the progress in the 
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understanding of the GRB phenomenon is enormous (e.g., Zhang & Meszaros 2004; Piran 
2005). From the observational point of view nearly every individual afterglow has turned 
out to be specific in some sense. The richness in afterglow properties observed so far is not 
much surprising, however, given the fact that the afterglow phenomenon combines internal 
properties of the underlying burster population with properties of the external circumburster 
medium. This in combination with redshift provides a large parameter space for various fla- 
vors of GRB afterglows. However, in spite of their individualities, as an ensemble afterglows 
trace the underlying physical boundary conditions and the parameter space of the physical 
processes involved. Revealing this parameter space is of fundamental interest since it might 
much improve the understanding of the afterglow phenomenon. This has already motivated 
several groups to perform a systematic analysis of observational data from various afterglows 
(e.g., Frontera et al. 2000; Panaitescu & Kumar 2001a,b, 2002; Yost et al. 2003; Frail et al. 
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2003a, 2004; Stratta et al. 2004; Gendre & Boer 2005; Panaitescu 2005a,b). The present 
work continues this kind of investigations but concentrates on the optical/NIR bands. 

While a fit of the observed broad-band spectral energy distribution (SED) of an afterglow 
from the radio to the X-ray band is the best way in order to extract the underlying physical 
parameters (e.g., Frail et al. 2003b; Yost et al. 2002; Panaitescu 2005b), one might expect that 
the optical/NIR data alone are completely explainable within the context of the underlying 
theoretical model. In particular, in most cases it is the optical/NIR region where the sampling 
and the quality of the afterglow data is best. It is clear that any reliable theoretical model 
must be able to explain the observed richness in the phenomenology of optical light curves 
of GRB afterglows. In this respect it is worth summarizing and exploring the available data 
base on afterglows in the optical/NIR bands. This holds in particular with respect to the 
new era in GRB research initiated by the launch and operation of the Swift satellite (Gehrels 
2004). 

In the present paper we continue our systematic analysis of afterglow parameters based 
on optical and near-infrared data (Zeh, Klose, & Hartmann 2004, hereafter Paper I; Zeh, 
Klose, & Hartmann 2005). Our study includes all afterglows in the pre- Swift era with 
sufficient published data. While in Paper I we analyzed the afterglow light curves with 
special emphasis on an underlying supernova component, the goal of the present paper is to 
explore the parameter space of the light curve parameters and of physical quantities related 
to them. In Paper III (Kann, Klose, & Zeh 2005), we expand this analysis to the spectral 
energy distribution of afterglows in the optical and near-infrared bands in order to search 
for signatures of dust in GRB host galaxies. 



2. Data analysis 

We collected from the literature all available photometric data on GRB afterglows ob- 
served by the end of 2004 and checked them for consistency. We modeled the light curve of 
the optical transient (OT) following a GRB as a composite of genuine afterglow light, super- 
nova (SN) light, and constant light from the underlying host galaxy. The afterglow is either 
modeled by a single or a smoothly broken double power-law according to Beuermann et al. 
(1999). Our fitting procedure is based on a x 2 minimization with a Levenberg-Marquardt 
iteration. This minimization technique provides formal uncertainties of the fitted parameters 
through the covariance matrix. For the remainder of this paper, all uncertainties are la. 

We always fitted photometric magnitudes with a fitting equation as follows: 
m OT (t) = -2.5 1og{10- a4mc [(t/t b ) Qin +(t/t b ) Q2n ]- 1 / n 
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The parameters in Equation (1) are the prebreak decay slope a±, the postbreak decay slope 
a2, the break time tb, the steepness of the break n, the brightness of the host galaxy mi, ost , the 
constant m c which corresponds to the magnitude of the fitted light curve for the case n = oo 
at the break time tb (it stands for the intersection point of the prebreak and the postbreak 
slope, without considering a smooth transition), as well as the supernova parameters k and 
s, which indicate the luminosity ratio and the stretch factor normalized to SN 1998bw (cf. 
Paper I). If there is no break in the light curve, then Equation (1) reduces to 



where mi is the brightness of the afterglow at t — 1 day after the burst. Before the fit- 
ting process the data were corrected for Galactic extinction using the maps of Schlegel, 
Finkbeiner, & Davis (1998). For a more detailed description of our procedure we refer to 
Paper I. 

In the present paper, our input list contains 59 optical/NIR afterglows observed in the 
pre- Swift era which allowed us to perform a fit (Table 1). Sixteen of these have a well- 
sampled light curve in at least one photometric band (mostly the Rc band). These light 
curves show a well-detectable break and have sufficient data points before and after the 
break time, so that the prebreak decay slope and the postbreak decay slope according to 
Equation (1) are well-defined. More precisely, our selection criterium for the best-defined 
afterglow light curves is that the la error is less than 0.2 for a 1; less than 0.3 for a 2 , and 
less than 1 day for t b (Table 1). We did not use the accuracy of the fit (x 2 per degrees of 
freedom) as a selection criterium, however. We consider Equation (1) as an empirical first 
order approximation of an observed light curve and all deviations from the corresponding 
fit as fine structure (§ 4.4). Light curves with no break are by definition excluded from this 
sample of best-defined afterglows, because it cannot be ascertained with complete certainty 
if the slope is a pre- or postbreak. Whenever we did not detect a break, the data quality was 
usually insufficient to exclude the possibility of a jet break in the light curve. In these cases 
there is either no early time data available, or the break could have been missed because of 
a bright host galaxy or an underlying SN component. 

All but one (GRB 980519) of the 16 afterglows in our sample have a known redshift, 
none of these bursts is classified as an X-ray flash 1 . For the present study the afterglows 



1 While GRB 030429 is classified as an XRF in the observer frame (Sakamoto ct al. 2005), it would be an 
X-ray rich GRB in the host frame with E pea ^ = 93^1 keV. 



mcyr(t) 



+ k 10"' 




(2) 



-4- 



of GRB 021004 and GRB 030329 were excluded from our analysis because of their many 
rebrightening episodes (see also § 4.4). In the case of the afterglow of GRB 030329, we 
performed two additional fits that are also given in Table 1, the details are in Appendix A. 
While these fits give very good results, they are based on only a small part of the light curve 
and we thus do not include their results in the selected sample for statistical study, with the 
exception of Figure 8. The afterglow of GRB 000301C also shows a rebrightening episode 
(around 3.5 days after the burst; Garnavich, Loeb, & Stanek 2000; Gaudi, Granot & Loeb 
2001), but as the deviations are not that large and occur only during a certain period, the 
light curve can still be fitted with a broken power-law, even though Xd.o.f. becomes relatively 
large. 

A special note is required for GRB 021211. The afterglow light curve of this burst 
fulfills the aforementioned criteria for the amount of the la error bars of the fit parameters, 
nevertheless we have not included it in our subsample of well-defined light curves. According 
to our data, the light curve shows a break 0.11±0.09 days after the burst with a steepening 
by Act = «2 — cui = 0.26 ± 0.11. This is in close agreement with the amount of steepening 
expected for the passage of the cooling break across the optical window (for a Compton 
parameter less than 1; e.g., Panaitescu & Kumar 2001a). In addition, Nysewander et al. 
(2005) report on evidence for color changes of the afterglow just around this time. Given the 
fact that no optical data have been reported in the literature for the time period between 1 
and 10 days after the burst, we consider it as very likely that the break in the light curve 
we have found does indeed signal the passage of the cooling break (as already suspected by 
Nysewander et al. 2005, based on their finding of color changes), while the real jet break 
occurred between 1 and 10 days after the burst. Note that the light curve break around 
0.11 days is not identical to the break discussed by others concerning the very early light 
curve, which has been attributed to the reverse shock (Li et al. 2003; Fox et al. 2003; Wei 
2003; Holland et al. 2004; Kumar & Panaitescu 2003). Given these findings, and since we 
try to keep our subsample as homogeneous as possible, we have not included this burst in 
the following study. 

The light curve parameters we have deduced from our data might be slightly different 
to those obtained and used by other groups. This is mainly due to the fact that we use a 
different, and most likely larger data base. In addition, there is still some bias in the selection 
of the data, in the definition of what are outliers, which data should be used and which data 
should not. While we do not claim that our data base is the best for every individual GRB, 
most probably it is the most comprehensive set. The strength of our approach is that we 
analyze all afterglow light curves using the same numerical procedure. In this respect we are 
confident that the results we obtain are statistically robust. 
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3. Results 

The results of our light curve fitting of the individual afterglows are summarized in 
Table 1. Here, the second column indicates the sample of the best-defined 16 afterglows. In 
some cases we could only derive upper or lower limits for some parameters, because the light 
curve was poorly sampled (see § A). In the following we discuss only the results obtained 
for this sample of 16 afterglows. 

The distributions of the deduced light curve parameters are shown in Figures 2 to 4. 
The mean of the prebreak decay slope is a\ = 1.05 ± 0.10, with ai ranging from 0.58 ± 0.05 
(GRB 000301C) to a x = 1.76 ± 0.05 (GRB 011121), while a 2 ranges from 1.30 ± 0.02 (GRB 
041006) to 3.03 ± 0.27 (GRB 030429) with the mean at 2.12 ± 0.14. About half of the 
postbreak decay slopes have a 2 < 2.0. Both distributions overlap in the interval from about 
1.3 to 1.7. The distribution of a 2 is almost constant, with a possible cutoff around 1.3 and 
no preference for any value. 

Figure 3 displays the distribution of the difference of the decay slopes before and after 
the break, Aa = ol 2 — ol\. The distribution is asymmetric with its maximum around 0.8 
and a longer tail towards higher values. It is notably broader than the distribution of a\ 
and «2- R also indicates the possibility that some afterglows could have very shallow breaks 
with Aa < 0.3 that could easily be missed. The gap around Aa = 1.7 is probably due to 
low-number statistics and we do not consider it to be significant. In this respect we cannot 
confirm the potential evidence for a bimodality of the distribution of Aa (Panaitescu 2005a), 
allthough we cannot reject this possibility either. 

In Figure 4 we present the distribution of the observed break times for the afterglows 
with the best-defined light curves, but translated into the corresponding GRB host frame. 
This distribution is strongly asymmetric with a clear peak in the host frame at lower values 
around 0.3 days. That most breaks occur at relatively early times supports the view that in 
several cases (Table 1) light curve breaks might have been missed due to a lack of early-time 
data. The afterglow with the earliest break in the host frame was GRB 041006, while the 
afterglow of GRB 000301C had the latest break time. On the other hand, late breaks might 
have been missed in several cases too, because the afterglow was already too faint at the 
break time, and/or an underlying SN component or a bright host galaxy simply made the 
discovery of the break in the optical bands impossible. We conclude that these data indicate, 
even though they do not prove, that in fact all afterglow light curves have detectable breaks 
due to a collimated explosion, as long as they are not hidden by rebrightening episodes as in 
GRB 030329 (e.g., Lipkin et al. 2004). It is clear that any model that explains the observed 
light curve breaks must be able to reproduce this observed distribution (Figure 4). 
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The distribution of the shape parameter n (Equation 1) is more difficult to quantify. 
While for about half of the afterglows the data allowed us to determine n during the fitting 
procedure, in the other cases n did not converge during the fit (n — > oo), because the 
sampling of the data is not good enough. It is worth noting that whenever we were able to 
determine n, we obtained a relatively soft break {n ~ 1 — 2) and in each case the prebreak 
decay slope cui was very shallow. In the other cases we had to fix n at a relatively large 
value in order to obtain an acceptable fit, and we chose n = 10. Choosing n « 1 — 2 made 
the fit worse. Whether this indicates a possible bimodal distribution of the parameter n is 
an open question. The distribution of the shape parameter n is the biggest unknown, so far, 
including its theoretical interpretation. 



4. Discussion 

4.1. Wind versus ISM models 

Figure 5 shows the correlation between the parameters a.\ and a 2 in comparison with 
eight standard afterglow models that cover the cases (1) v < is c , (2) v c < is, Y < 1, (3) 
v c < v, Y > 1, 2 < p < 3, (4) v c < v, Y > l,p > 3 for the ISM and for the wind model, where 
v is the observed frequency and Y is the Compton parameter (Panaitescu & Kumar 2001a, 
their eqs. 21, 22 and 29). Those models with v < v c and Y < 1 require p < 2 in order to 
explain an observed a 2 < 2. In particular, there is a group of five bursts (GRBs 990123, 
991216, 010222, 030328, 041006) that cluster around a 2 = 1.5. Within the corresponding la 
error bars a 2 > 2 is basically ruled out and so is p > 2. Such shallow postbreak decay slopes 
cannot be explained by a flat electron distribution either (Dai & Cheng 2001; Bhattacharya 
2001; Wu et al. 2004). 

Based on the underlying theoretical models, which predict p = a 2 , then we find that 
the parameter space of p is rather broad, ranging from about 1.5 to 3. While the results 
obtained for p for the individual afterglows differ among various authors, all studies agree 
that within the current theoretical framework no evidence for a universality of p is found 
(e.g., Panaitescu & Kumar 2001a,b, 2002; Preece et al. 2002; Yost et al. 2003; Panaitescu 
2005b, ;Paper III). This is contrary to what one might expect from theoretical models of 
highly relativistic shocks (Achterberg et al. 2001; Kirk et al. 2000), and contrary to what 
one might prefer on theoretical grounds (Freedman & Waxman 2001). 

If one allows for p < 2 then an inspection of Figure 5 shows that a wind model with 
v < u c is preferred for GRBs 980519, 990123, 991216, 000926, 020124, 020405, 030328, and 
041006, even though GRBs 990123, 991216 and 000926 lie fairly off the models because of 
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their relatively small a 2 or large a\. On the other hand, an ISM model with v < v c is 
preferred for GRB 990510 and GRB 011211. GRB 020813 is either an ISM case with v < v c 
or an ISM/wind case with v > u c , in the following, we will assume this to be an ISM case. 
Two afterglows (GRB 010222 and GRB 011121) could be an ISM or a wind case with v > v c . 
In the case of GRB 011121, the error bars are so large that v < v c cannot be fully excluded. 
Finally, the afterglows of GRBs 000301C, 030226, and 030429 are outliers because of their 
relatively large o 2 - It is noteworthy that all afterglows with soft breaks in the light curve 
(n ~ 1 — 2) belong to the group of bursts that are less compatible with a wind profile, and 
all these afterglows have a± < 1. 

It is difficult to quantify whether the outliers really represent a different population or 
whether in these cases the light curves are simply ill-defined. For example, the afterglow of 
GRB 00030 1C was affected by a strong rebrightening episode that has been modeled by a 
gravitational microlensing event (Garnavich, Loeb, & Stanek 2000). As these authors note, 
the removal of this event leads to «i ~ 1.1, which would shift GRB 00030 1C toward the 
theoretical prediction of the ISM model with v < v c . In the case of GRB 030226 there 
is evidence that the afterglow light curve showed fluctuations. In combination with the 
relatively sparse set of postbreak data it is well possible that the late-time observations 
stopped when the afterglow underwent a fluctuation, so that finally the deduced postbreak 
decay slope is too large. On the other hand, early-time spectra of this afterglow reveal 
features that can be best understood as due to a stellar wind profile (Klose et al. 2004). In 
the case of 030429, Jakobsson et al. (2004b) find that the light curve undergoes a significant 
rebrightening around 1.7 days after the burst. They suggest to exclude the data around the 
rebrightening from the fit and fix «2 = 1.7 (deduced from the SED and the a — (3 relations). 
In this case the light curve would be compatible with an ISM/wind model and v > u c , 
although it is unclear how large the error in a 2 would be. 

If we neglect GRBs 000301C, 030226, and 030429, then Figure 5 shows that the group 
of optical afterglows that is compatible with a wind model is notably larger than the group 
of afterglows that prefers an ISM model. While basically all studies in the literature agree 
that afterglows seem to separate into a group that can be best described by a wind model 
and a group that can be best described by an ISM model, our data show that the wind 
scenario is statistically preferred. In fact, within their la error bars nearly all (0:1,0:2) pairs 
are compatible with a wind profile while for the ISM model such a statement is clearly ruled 
out. 
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4.2. The jet opening angles 

Figure 6 displays the distribution of the jet half-opening angle, #j Ct , for our sample of 
16 bursts (GRB 980519 is not included as no redshift of this burst is known) as derived from 
the observed break time, assuming the uniform jet model (Rhoads 1999). We calculated #j et 
following Sari, Piran, & Halpern (1999) for an ISM medium (GRBs 990510, 000301C, 011211, 
020813, 030226, 030429) and Bloom, Frail & Kulkarni (2003) for a wind-like medium (GRBs 
990123, 991216, 000926, 010222, 011121, 020124, 020405, 030328, 041006), according to the 
results obtained for the density profile of the individual afterglows (§ 4.1). For the isotropic 
equivalent energy E iso , the radiation efficiency and the redshift we adopted the values given 
by Friedman & Bloom (2005). In the case of an ISM model we used the circumburster density 
as given in Friedman & Bloom (2005), while for the wind model we assumed a mass-loss 
rate to wind speed ratio of A* = 1 (cf. Chevalier & Li 2000). The distribution of 9j et that we 
have found is strongly asymmetric with a peak between 2 and 5 degrees, has a lower cutoff 
around 2 degrees and rapidly falls towards larger angles, in agreement with what has been 
found in previous studies (Frail et al. 2001; Panaitescu & Kumar 2001b, 2002; Bloom, Frail 
& Kulkarni 2003). 



4.3. Correlations 

Using the derived jet half-opening angles (§ 4.2) we find that the distribution of the 
beaming corrected energy release in the gamma-ray band ranges from E cor = 10 49 9 erg 
(GRB 041006) to E COI = 10 514 erg (GRB 990123). In combination with the corresponding 
peak energies, -E pea k, in the gamma-ray band (Friedman & Bloom 2005), in Figure 7 we 
plot the correlation between E cor and E pea k in the GRB host frame, as it was first reported 
by Ghirlanda, Ghisellini & Lazzati (2004). Considering bursts #9 and #12 as outliers and 
excluding them from the fit, we find -E pea k — 748 (E 5 i tCor ) v , with rj = 0.79 ± 0.09. This is 
in qualitative agreement with Ghirlanda et al. (2004) as well as Dai, Liang & Xu (2004). 
On the other hand, there are differences, in particular concerning the existence of the two 
outliers. They can be understood, however, since the fit includes assumptions about the 
gas density in the GRB environment, which enters the calculation of the jet opening angle. 
We made use of the values provided by Friedman & Bloom (2005) and these are higher by 
a factor of 10/3 compared to the values used by Ghirlanda et al. (2004). In addition, in 
several cases the gamma-ray data given in Friedman & Bloom (2005) are notably different 
from those used by Ghirlanda et al. (2004). Additionally, we also considered a wind-like 
circumburst medium for some cases to calculate the jet half-opening angle, while Ghirlanda 
et al. (2004) only regard an ISM-like circumburst medium. It is therefore not surprising 
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that we do not exactly reproduce their results (for a discussion, see also Friedman & Bloom 
2005). For instance, if we reduce the assumed circumburster gas density for bursts #9 and 
#12 from 10 to 1 cm -3 , then the corresponding data points do not fall out of the sample 
anymore. A fit then provides r) = 0.78 ± 0.09. While our data base is too small in order 
to investigate the role of outliers in the Ghirlanda relation, we note that also the relation 
between the isotropic equivalent energy release in the gamma-ray band and the intrinsic peak 
energy (Amati et al. 2002) is not very tight. In a recent study on BATSE GRBs, Nakar & 
Piran (2005) find that about 25% of all bursts do substantially deviate from this empirical 
relation (see also Band & Preece 2005). 

In addition to the Ghirlanda relation, we have searched for linear correlations between all 
individual afterglow parameters and between the burst parameters in the gamma-ray band 
and the corresponding afterglow parameters. Table 2 lists the corresponding correlation 
coefficients derived from weighted linear fits. The Ghirlanda relation (here including the two 
outliers discussed above) is between log(i? cor ) and \og(E pcak ). Relatively tight correlations 
between 9j et and h/(l + z) and between 6 ]et and log-E cor are expected, as 6- ]Ct derives from 
tb/(l + z) and log-E C or derives from #j Ct . Note that the correlations for the wind model are 
much tighter than for the ISM model, giving further significance to the statistical conclusion 
that most circumburst environments are wind-blown (§ 4.1). Next to the Ghirlanda relation 
and the others just discussed, there are more correlations that seem significant according to 
the absolute correlation coefficient, but visual inspection does not support this conclusion. 
This holds also for the potential correlation between 9- jCt and a x which has been reported for 
X-ray data (Liang 2004). Still, it is interesting to note that a 2 seems more or less correlated 
with all other parameters, including the redshift. On the other hand, 0j et , log-E^or and 
log-Epeak are completely uncorrelated with the redshift. We conclude that no evolutionary 
effect in the initial explosion parameters is evident over a wide range of redshifts. 

Of special interest is the parameter n, which indicates the smoothness of the break 
(Equation 1). Even though we could only determine n for a few afterglows, it looks suspicious 
that each time we had to fix n to a relatively high value in order to get an acceptable fit, 
the prebreak decay slope is a\ > 1. To test if this could be due to a numerical problem we 
reconsidered the afterglow light curve of GRB 030329 and fitted it only between 0.28 and 1 
days, which includes the time around the supposed jet break (e.g., Uemura et al. 2003) but 
excludes the cooling break (Sato et al. 2003) and the rebrightening episodes (e.g., Lipkin et 
al. 2004). For this time period the data density is high enough in order to deduce a value for 
n even if the break is very sharp. Indeed, in this case we find «i = 1.17±0.01, t b = 0.68±0.02 
days, a 2 = 2.21 ± 0.07 and n = 7.54 ± 1.47. Adding this to our sample of afterglows with a 
deduced parameter n (Table 1), a weak trend between n and a\ becomes apparent (Fig. 8). 
It indicates that a shallow prebreak decline leads to a smooth break or, seen the other way 
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around, a smooth break (a "rollover") implies a shallow prebreak decay slope. As we only 
have very few values of n, we do not regard this as strong statistical evidence for a correlation, 
since an observational bias cannot be excluded. Afterglows with a shallow decay are bright 
for a longer period of time, which makes them easier to follow. Therefore, the data density 
around the break time is usually higher compared to most afterglows with a steep decline. 
A high data density around the break time is essential to determine a value for n, however. 
Nevertheless, it is worth to check if this trend is confirmed in the Swift-era. 

Unfortunately, most bursts in our sample are at high redshift, so that no supernova data 
are available. Consequently, no statistically founded conclusions can be drawn on a poten- 
tial correlation between afterglow parameters and the corresponding supernova parameters 
(Paper I). 

4.4. Fine structure in the light curves 

In principle, Equation (1) represents a first order approximation of an observed afterglow 
light curve. We consider any residuals that remain after subtraction of the fit from the 
observational data as the fine structure in the light curve. Since the detection of fine structure 
in optical afterglows depends strongly on the sampling density of the light curve and the 
quality of the data, an observational bias cannot be excluded, which makes it difficult to 
compare the fine structure among the individual afterglows. On the other hand, it is well- 
known that some afterglows show basically little or no evidence for fine structure when 
sampled very densely with the same telescope (the best example being GRB 020813; Laursen 
& Stanek 2003), while others show a substantial amount of fine structure (e.g., GRB 030329; 
Lipkin et al. 2004, their figure 4). In order to investigate the occurrence of fine structure in 
more detail, we have shifted all residuals to a common evolutionary phase of an afterglow. 
We favor the idea that this can be done by normalizing the time t that has elapsed sind the 
GRB trigger to the break time t b of the corresponding afterglow (Table 1). While Figure 1 
displays the individual light curves and corresponding residuals that remain after subtraction 
of the fitted curves from the observational data of the 16 afterglows in our sample, Figure 9 
displays all residuals in a single plot as a function of t/tb- We have included here only those 
data with less or equal 0.05 mag individual photometric error. The ratio t/tb is independent 
of redshift and allows us to draw some general conclusions about the occurrence of fine 
structure in afterglow light curves. 

First at all, again we find no evolutionary effect in the data. The width of the magnitude 
distribution of the fine structure of all 16 afterglows in the prebreak evolutionary era spanning 
one decade in time (0.1 < t/t b < 1), is identical to the width of the magnitude distribution 
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in the postbreak era spanning one decade in time (1 < t/4 < 10), namely ±0.1 magnitudes. 
The handful of data points around t/t}, = 0.2 that reach beyond ±0.2 magnitudes mainly 
belong to GRB 011211 (cf. Holland et al. 2002; Jakobsson et al. 2004a) and is statistically 
not significant. This picture does not change if we allow for larger individual photometric 
errors but then the statistical significance of this finding becomes less strong. We conclude 
that, on average, a patchy surface structure of afterglow shock fronts (Meszaros, Rees, & 
Wijers 1998; Nakar, Piran & Granot 2003) is not present in the photometry at times later 
than t/t b > 0.1, with the probable exception of GRB 011211. In addition, we find that the 
amplitude of the fine structure of all 16 afterglows as a group (±0.1 mag) is smaller by a 
factor of 4 compared to the fine structure (or fluctuations) seen in the optical afterglows of 
GRB 021004 (e.g., de Ugarte Postigo et al. 2005) and GRB 030329 (Lipkin et al. 2004), which 
are plotted in comparison 2 . In other words, the latter two optical afterglows are indeed very 
different from all other afterglows we have investigated, as the deviations persist even after 
the break. 

On the other hand, the residuals of the early afterglows of GRB 021004 and GRB 030329 
are very similar. We find via x 2 minimization that a shift of the GRB 021004 light curve in 
t/h by a factor of about 2.7 superposes the early light curve evolution of both bursts (the x 2 
minimum we find is not sharp, shifts between 2.4 and 3.0 are acceptable). This is astounding, 
as the bursts happened at two very different cosmological epochs. A deeper analysis of this 
result will be pursued in a future publication (Kann et al. 2005, in preparation). 

5. Summary 

Based on a systematic analysis of the optical and NIR data of all GRB afterglows with 
sufficient published data in the pre-Swift era we have explored the parameter space of the 
afterglow light curves and of physical quantities related to them. From the 59 afterglows 
investigated (Table 1) we constructed a sample of 16 bursts with the best-defined light curves 
useful for our purposes. Thereby we excluded the afterglows of GRB 021004 and GRB 030329 
because of their many re-brightening episodes which made it difficult to fit them. 

Using the sample of the 16 afterglows with the best defined light curves, we find that 
in the optical bands, the average afterglow light curve is characterized by a prebreak decay 
slope ai = 1.0 ± 0.1 and a postbreak decay slope a 2 = 2.1 ± 0.1. The distribution of both 
parameters is rather broad but possible cutoffs are apparent in the available data set. In 



2 Note that our broken power law fits find late breaks for both these afterglows (Table 1), so most of the 
data is at t/tb < 1- 
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particular, there is no evidence for a universality of a 2 as it has been predicted in some 
afterglow models. The distribution of the break time in the host frame rises sharply towards 
smaller values, with the most likely value at t b / (1 + z) — 0.3 ± 0.2 days. 

We have then used the deduced light curve parameters to extract information about 
the nature of the GRB environment using standard afterglow models (Panaitescu & Kumar 
2001a). We find that in most if not all cases the data are in agreement with a wind model. 
A general preference of an ISM model is ruled out. In addition, we find that the distribution 
of the power-law index p of the electron distribution function is rather broad ranging from 
about 1.5 to 3, supporting the view of a non-universality of p, in agreement with other studies 
(e.g., Panaitescu & Kumar 2001a,b, 2002; Preece et al. 2002; Yost et al. 2003; Panaitescu 
2005b, Paper III). Furthermore, we have searched in our data set for potential correlations 
between the various light curve parameters and those that characterize the corresponding 
burst in the gamma-ray band. With the exception of the Ghirlanda relation between the 
beaming corrected energy release in the gamma-ray band and peak energy in the GRB host 
frame (Ghirlanda et al. 2004), no other tight correlation has been found. An intriguing 
correlation may exist between the pre-break decay slope a\ and the smoothness parameter 
of the break n, but more data is needed to verify this. 

Finally we have analyzed in which manner the data indicate to a general fine structure 
that is superimposed a light curve decay according to the empirical Beuermann double power- 
law. When normalized to the corresponding break time t b of a burst we find no evidence 
in the data that there is more structure in the light curves at times 0.1 <£/£&< 1 than 
at times 1 < t/t b < 10. On the other hand, we find that the afterglows of GRB 021004 
and GRB 030329 are very different from all 16 afterglows in our sample. While the latter 
vary on average by only 0.1 magnitudes around the fitted light curve, the former vary by 
0.4 magnitudes. Moreover, the fine structure of the light curves of GRB 021004 and GRB 
030329 are initially very similar. 

It is clear that more afterglows with well-sampled optical light curves are needed in 
order to get better insight onto the parameter space of the physical processes involved. Even 
in the Swift era this might not be an easy task since most afterglows are simply very faint 
some days after the burst. The observational challenge therefore remains the availability of 
observing time on large optical telescopes in order to determine the light curve parameters 
as well as possible. 
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S.K. acknowledge financial support by DFG grant Kl 766/13-2. We wish to thank S. Covino, 
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thank the anonymous referee for very constructive remarks. 

A. Notes on individual bursts 

GRB 970508 The afterglow light curve of this GRB is anomalous, featuring an early 
plateau phase enduring until one day after, followed by a steep (a ~ —3.4) re-brightening. 
Starting at 1.9 days, the afterglow decays with a simple power law. Our fit starts at this 
point. The light curve is well sampled and no break is seen, thus it is possible that the decay 
is postbreak (the break being hidden by the early anomalous behavior). 

GRB 970815, GRB 030131 In each case, the afterglow light curve has only two data 
points, and a late upper limit indicates a faint host. As the degree of freedom is zero, no 
errors are given. 

GRB 980326 The light curve fit includes a supernova component with k — 1, s — 1 fixed 
(equation 1), assuming a redshift of z — 1 (Bloom et al. 1999). 

GRB 980519 The first data point of the Rc band light curve of the afterglow of GRB 
980519 is at t ~ 0.5 days after the burst. In the Ic band and V band earlier data points 
exist but no late-time data. In order to improve the fit, we assumed achromacy and mixed 
these bands by shifting the Ic and V band to the same zero point as the Rc band, fitting 
the composite light curve. For + z), we assumed a redshift of z — 1.5, since Jaunsen 
et al. (2001) state that z > 1.5 from the absence of a supernova bump in the light curve. 
Using z = 1.5 also gives a very good SED fit (Paper III). 

GRB 990123 and GRB 021211 These afterglow light curves have very early detections, 
where the light curve is dominated by reverse shock emission. Data from the reverse shock 
dominated phase have not been included in the fits. For GRB 021211 the result are sensitive 
to the data used for the fit. We used only data with t > 0.014 days. See also our comment 
in §2. 
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GRB 990705, GRB 020322, GRB 020410, GRB 030115, GRB 030324, GRB 
031220, GRB 040422 In all these cases, afterglow data are too sparse to confine certain 
parameters of the light curves. The addition of observational upper limits was however used 
in order to derive upper or lower limits on these parameters. For light curve decay slopes, the 
sequence data point - upper limit - host magnitude sets a lower limit on the decay rate (e.g., 
GRB 020410) as the observational upper limit does not preclude an even steeper decay rate. 
On the other hand, the sequence upper limit - data point - host magnitude sets an upper 
limit on the decay rate (e.g., GRB 030324), as the upper limit cannot preclude a slower 
decay rate. In the case of GRB 040422, early unfiltered upper limits were, after correction 
for Galactic extinction, brightened by three magnitudes (the approximate typical Rc — K 
color) to get early .fT-band upper limits. In some cases early data was adequate to derive ot\ 
and a later upper limit lies beneath the extrapolated light curve decay, indicating a break 
in the light curve must have occurred. 

GRB 991208, GRB 000131, GRB 001007 In all three cases, the optical afterglows 
were located several days after the burst and exhibited a steep decay a ^ 2. It is highly 
probable that the first observations are after the jet break, meaning that the decay slope is 
«2. For GRB 991208, the upper limit on a\ stems from a very shallow but very early upper 
limit (Castro-Tirado et al. 2001). We note that a postbreak decay slope is also possible for 
GRB 000911 (Masetti et al. 2005) - if so, it would be quite shallow. There are also several 
GRBs (e.g., GRB 990308, GRB 001011) where data quality is so sparse that no evidence for 
or against a break in the afterglow light curve can be found. 

GRB 011121, GRB 020405 These bursts have relatively bright or structured hosts, 
which makes it difficult to extract the afterglow light curve. Different groups used different 
methods to do this. Therefore we decided not to mix the late data points (where the different 
host galaxy subtraction methods lead to significantly different magnitudes), but instead to 
use the data set from only one group in each case. The late data set of GRB 011121 is 
extracted from Greiner et al. (2003), the data set from GRB 020405 is provided by N. 
Masetti (private communication). 

GRB 020305, GRB 030725 The late-time rebrightenings of these afterglow light curves, 
which have been attributed to possible supernova components (Gorosabel et al. 2005a; 
Pugliese et al. 2005), have not been included in the fit as no reliable redshift estimate is 
known. In the case of GRB 030725, there is a large data gap from 0.5 to 4.5 days. A fit that 
leaves t b as a free parameter is not confined in oi\ and t b , thus, t b was fixed at a reasonable 
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value. 

GRB 021004, GRB 030329 The light curves of these afterglows are clear outliers com- 
pared to the rest of our sample and very complicated to fit. The light curves each show 
several re-brightening episodes, which cannot be fit correctly with a smoothly broken power- 
law. As especially the light curve of the afterglow of GRB 030329 has been analyzed in 
detail (e.g., Lipkin et al. 2004), we did not do this again. As with the other afterglows, we 
fit smoothly broken power-laws to the light curves (the results are given in Table 1). In the 
case of GRB 021004, the anomalous behavior at t < 0.07 days is excluded from the fit. We 
performed these fits for completeness, even though the Xd.o.f. values (xio.f. > 20 for both) 
show that they are not well approximated by equation (1). While a± is almost unaffected, 
a 2 and tb are highly dependent upon the value of the smoothness parameter n, which had 
to be fixed. These results are thus not included in our statistical analysis. 

If we concentrate on certain parts of the light curve of the afterglow of GRB 030329, they 
can be very well fit by a smoothly broken power-law. We have performed two additional fits, 
one using only data up to 0.55 days, before the time of the probable jet break (e.g., Uemura 
et al. 2003), and thus encompassing the cooling break (Sato et al. 2003), the other one using 
data from 0.28 days up to 1 day after the burst (the beginning of the first rebrightening) 
and thus encompassing the supposed jet break. In the former case, n is not confined, the 
fit formally finds n = 450 ± 2800, but a very sharp transition is indicated, and thus we fix 
n = 100. In the last case, we were able to let the smoothness parameter n vary and obtain 
n = 7.54 ± 1.47 for the jet break. The values we derive for the jet break, a\ = 1.17 ± 0.01 
and a 2 = 2.21 ± 0.07, are both unremarkable and lie well within the distribution we find for 
the 16 afterglows we study here (cf. Figure 2). The slope change Aa is Aa = 0.33 ±0.01 for 
the cooling break (slightly higher than the theoretical prediction of Aa = 0.25, Panaitescu 
& Kumar 2001a) and Aa = 1.04 ± 0.10 for the jet break. This value is also unremarkable 
(cf. Figure 3). The rest frame break time is 0.58 ± 0.01 days after the burst, once again 
typical for our sample of 16 afterglows (cf. Figure 4). The values of a± and a 2 lie very well 
on the theoretical line describing a Wind/ISM model with v > v c . It is one of only a few 
bursts that are found in this region (cf. § 4.1, Figure 5). As the cooling frequency has passed 
the optical bands before the jet break and now lies at longer wavelengths, this implies an 
evolution from high to low frequencies and thus a wind model (Chevalier & Li 2000). 

GRB 030323 The afterglow of this GRB has a late break which is only represented by 
one data point. a 2 was fixed to a value derived from a free fit to derive meaningful errors 
on a\ and 
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XRF 030723 The early light curve of this afterglow has a plateau phase. We do not fit this 
plateau phase but start at 0.9 days after the burst. While a light curve break is evident, the 
prebreak data are inadequate to give more than limits. We also do not fit the late time data 
points (t > 8 days), which show a significant rebrightening (Fynbo et al. 2004) which can be 
modelled with a SN light curve only when it rises much steeper than SN1998bw. Another 
explanation could be that this rebrightening is caused by a two component jet (Huang et al. 
2004). 

GRB 040827 The steep decay of this light curve (a > 2) indicates that this is a postbreak 
decay, but the data quality is low and the error large enough to make this conclusion unsure. 

GRB 040924 The light curve is fit with an unbroken power-law, excluding the earliest 
data point (Fox 2005) from the fit. 

GRB 041006 As the photometric calibration of the very earliest data point of the optical 
afterglow (Maeno et al. 2004) is unsure, it is not included in the fit. Including it would 
strongly reduce «i and n. 



B. References of host magnitudes 

The following magnitudes have not been corrected for Galactic extinction and thus differ 
from those displayed in Table I. 

• GRB 990308: R c = 29 A ± 0.4 mag, Jaunsen et al. (2003) 

• GRB 001011: R c = 25.38 ± 0.25 mag, Gorosabel et al. (2002) 
GRB 020305: R c = 25.17 ±0.14 mag, Gorosabel et al. (2005a) 

• GRB 021004: R c = 24.21 ± 0.04 mag, de Ugarte Postigo et al. (2005) 

• GRB 030115: R c = 25.2 ± 0.3 mag, Dullighan et al. (2004) 

• GRB 030324: R c = 25.16 ± 0.24 mag, Nysewander et al. (2004) 

• GRB 030329: R c = 22.66 ± 0.04 mag, Gorosabel et al. (2005b) 
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Table 1. Results of the light curve fitting (Equation 1) 
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Table 1 — Continued 
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Table 1 — Continued 
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a If a host magnitude is given with an error then this is the result of the fit. Otherwise we had to fix this value 
because the data set at late times is too sparse. In such cases we either used the host magnitudes reported in 
the literature (GRBs 990308, 001011, 020305, 021004, 030115, 030324, 030329; see § B) or we used a reasonable 
estimate (GRBs 970815, 990705, 000131, 030131, 030226, 0303227, 030418, 030429, 030723, 030725, 040106, 
040916, 041006). 



b See § A. 

c This fit uses only data up to 0.55 days after the burst, and encompasses the probable cooling break (Sato et 
al. 2003). See § A. 

d This fit uses only data from 0.28 days to 1 day after the burst, and encompasses the supposed jet break (e.g., 
Uemura et al. 2003). See § A. 



Table 2. The absolute value of the correlation coefficient between the various parameters 
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a Here the data are calculated according to an ISM or a wind model depending on the location of the corresponding burst in 
the cki — a 2 -diagram (Fig. 5 and § 4.2). 
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Fig. 1. — The light curves of the 16 afterglows we have investigated. In all but one case (GRB 
990510) we have analyzed the ifo-band data. Am is the difference between the observed and 
the fitted magnitude. For GRBs 011121, 011211, and 020405 the data are corrected for the 
flux from the underlying host galaxy. For GRBs 030226 and 041006 the -R^-band magnitude 
of the host was assumed to be 27.1 and 28.4, respectively. The numbering follows Table 1. 
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Fig. 1. 



continued. 
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Fig. 2. — The distribution of the prebreak decay slope «i and the postbreak decay slope 
a 2 for the 16 afterglows in our sample with the best defined light curves (Table 1; bin size 
0.5). The broken line and the solid line represent the cumulative distribution for a.\ and a 2 , 
respectively. 




Fig. 3. — The distribution of the light curve steepening Aa = 0,2 — ol\. The numbering 
follows Table I. The solid line is the cumulative distribution. 
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Fig. 4. — The distribution of the break time tb in the GRB host frame for all bursts of our 
sample of 16 afterglows. The data cover the range from %/ (1 + z) = 0.14 ± 0.02 days (GRB 
041006) to 1.62 ± 0.06 days (GRB 000301C). Most afterglows exhibit a break at less than 1 
day after the burst in the host frame. For GRB 980519 we assumed a redshift of z — 1.5, 
but it would fall into the first bin even if z — 0. The solid line represents the cumulative 
distribution. 
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Fig. 5. — The afterglow parameters ot\ and a 2 i n comparison with standard afterglow models 
(Panaitescu & Kumar 2001a). The dashed line is for an ISM density profile with v < v c , 
the dotted line for the ISM/wind model with v > v c and a Compton parameter, Y, less 
than 1, and the solid line for the wind model with v < v c . Note that we have extended 
these curves to a 2 < 2.0. The two curves around (cti = 1.0, at 2 = 2.0) represent the case 
v > v c , Y > 1,2 < p < 3 (Panaitescu & Kumar 2001a). The dash-dotted line stands for 
the theoretical prediction of the passage of the cooling break (a 2 — cxi + 0.25). While GRB 
990123 (# 2) lies on this line within errors, it has a late achromatic break which is very 
probably a shallow jet break. 
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Fig. 6. — The distribution of the derived jet half-opening angle, 6>j Ct , of our sample (§ 4.2). 
GRB 980519 is not included here since its redshift is not exactly known. 
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Fig. 7. — The Ghirlanda relation (Ghirlanda et al. 2004) as it follows from our light curve 
data (Table 1, § 4.2) in combination with the high-energy data given in Friedman & Bloom 
(2005). GRB #1 is not shown because of the uncertainty of its redshift and GRBs #5, #6, 
jfl are not included here since E ves ^ is not known. Our fit gives a slope of 0.79 ± 0.09 if the 
outliers GRB 011211 (# 9) and 020813 (# 12) are excluded. For comparison, the dotted 
line shows the relation obtained by Ghirlanda et al. (2004) based on their data base. 
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Fig. 8. — The relation between the prebreak decay slope ol\ and the smoothness of the break 
n (Equation 1). For all afterglows not in this figure we had to fix the smoothness parameter 
to a relatively high value of n — 10 to get an acceptable fit (Table 1). All these fits also 
find a-i > 1. In addition, we include here a special fit of the afterglow of GRB 030329 as 
explained in the text (§ 4.3 and Appendix A). This is the only fit where a\ > 1 and a value 
for n could be deduced. 
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Fig. 9. — The residuals Am = m b seive d — rngt for all 16 afterglows in our sample (Table 1, the 
large black points). Included are only those data points where the reported magnitude error 
is less or equal than 0.05 mag. The time axis is normalized to the corresponding break time 
of the burst. The data indicate that any fine structure in the light curves is on average less 
than ±0.1 magnitudes with no evidence for evolution. Note that the ratio t/tj, is independent 
of redshift. For comparison, the residuals of GRB 021004 (upward pointing blue triangles) 
and GRB 030329 (downward pointing red triangles) are plotted. The systematic deviations 
from the Beuermann law are clearly seen, reaching almost 0.4 magnitudes. Furthermore, 
the initial fine structure of GRB 021004 and GRB 030329 is very similar; shifting the GRB 
021004 light curve by a factor of 2.7 superposes this initial fine structure. 



